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1 .  Introduction.  If  a  random  variable  is  exponentially  distributed  with 

2  2 

y  *  EX  and  y_  *  EX  ,  then  y,  *  2y  .  Defining  p  =  [ — =--l| ,/  it  is 

i  2yZ  • 

tempting  to  conjecture  that  under  mild  restrictions  a  distibution  with 

small  p  is  approximately  exponential.  That  restrictions  are  needed  is 

1  - ) 

seen  by  the  example,  Pr(X»0)  »  Pr(X-l)  *  for  which  p  =  0. 

The  scale  invariant  quantity,  p,  was  suggested  by  Keilson  [7]^.  It 

—2.  x  m 

has  an  interesting  interpretation.  Define  G(x)  *  y  F(s)ds,  the 

stationary^  Renewal  distribution  corresponding  to  F.  Then  uG  *  m2/2p 
UG 

and  p  a|  - -  l|.  The  parameter  p  is  thus  the  scaled  (by  u)  distance 

between  y  and  y^.  For  F  exponential,  F  *  G  and  thus  y  ■  Pg- , 

The  problem  of  interest  can  be  stated  as  follows:  Given  a  class 

X  of  distributions,  along  with  the  first  two  moments  y  and  y^  find 

upper  bounds  for  sup  {F(t)-e  in  terms  of  p. 

v  t.Fei 

The  above  problem  for  the  class  of  completely  monotone  distributions 
(mixtures  of  exponential  distributions)  was  studied  by  Keilson  [7],  Heyde 
[5],  Heyde  and  Leslie  [6],  Hall  [4],  and  Brown  [1], 

Brown  [1]  considered  the  class  of  IMRL  (increasing  mean  residual  life) 
distributions  on  [O,00)  deriving: 


(1.1) 

(1.2) 

(1.3) 

(1.4) 


sup|F(t)-e_t/yj  <^- 

sup | F(B)-G(B) j  <  -£r 
Beg  P 

sup|G(B)  -  y_1e"t/wdt:j  <  --r  \ 

Beg  .  JB  ^ 

-t/vr  n 

sup|G(t)-e  |  • 
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In  (1.2)  and  (1.3)  above  8  is  the  collection  of  Borel  subsets  of 

[  0 ,°°) .  The  quantity  was  shown  to  be  the  best  upper  bound  for  (1.1) 

p+l 

and  (1.2)  even  within  the  subclass  of  completely  monotone  distributions. 

Brown  [2]  considered  the  class  of  IFR  (increasing  failure  rate) 
distributions.  It  turns  out  that  in  this  case  (1.1)  and  (1.2)  hold  with 


replaced  by  2p,  and  (1.3)  and  (1.4)  hold  with 


replaced  by 


p.  The  bound  2p  is  the  best  bound  for  (1.1),  among  bounds  of  the 
form  cpa. 

In  this  paper  we  consider  the  problem  for  F  NBUE  (new  better  than 
used  in  expectation)  and  for  F  NWUE  (new  worse  than  used  in  expectation) . 
These  are  the  weakest  among  the  commonly  studied  classes  of  aging  distribu¬ 
tions,  and  it  is  often  easy  to  demonstrate  that  a  distribution  belongs  to 
one  of  these  classes  (NBUE  and  NWUE  are  defined  in  Section  2).  The  methods 
of  Brown  ([1],  [2])  do  not  generalize  to  these  cases  because  the  partial 
ordering  between  F  and  G  is  too  weak.  Instead  we  use  Fourier  methods 
adopted  from  Feller  [3].  Our  main  result  is  that  for  F  NBUE  or  NWUE: 


(1.5) 


sup|F(t)-e  <  Ap*^ 


where  A 


3.119.  For  the  NBUE  case  we  show  that  the  best  bound  of 


the  form  cpa  has  a  *  1/2  and  1  £  c  £  .  Thus  the  potential 

improvement  in  (1.5)  for  F  NBUE  is  the  lowering  of  the  constant  from 
3.119  to  1.  This  remains  true  even  within  the  subclass  of  IFRA  distri¬ 


butions. 
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2.  Definitions  and  Preliminary  Results.  A  distribution  F  on  [0,“) 
with  F(0)  <  1  and  finite  mean  y  is  defined  to  be  NBUE  if 
E(X-t | X>t)  <  u  for  all  t  >  0  with  F(t)  >  0.  Since  E(X-tjX>t)  * 
yG(t)/F(t),  it  follows  that  F  is  NBUE  if  and  only  if  F  is  stochas¬ 
tically  larger  than  G,  the  stationary  renewal  distribution  corresponding 
to  F.  Define  hg  to  be  the  failure  rate  function  of  G  and  note  that 
hg(t)  ■  [E(X-tjX>t>]  \  thus  F  is  NBUE  if  and  only  if  hg(t)  £  y  ^ 
for  all  t  ^  0  with  F(t)  >  0. 

A  distribution  F  on  [ 0 , °°)  with  F(0)  >  1  and  finite  mean 
is  defined  to  be  NWUE  if  E(X-t|x>t)  >  y  for  all  t  £  0  with  F(t)  >0. 

This  is  equivalent  to  F  being  stochastically  smaller  than  G,  and  also 

to  hG  <  y  1. 

Lemma  2.1.  If  F  is  NBUE  then  G(t)  <_  e  for  all  t  _>  0;  for  F  NWUE, 

G(t)  2.  e~t/u  for  all  t  >  0. 

Proof,  For  F  NBUE  let  tg  be  the  smallest  number  such  that  F(tg)*0, 

with  tQ  ■  00  if  F(t)  >  0  for  all  t.  Now  hg(t)  £  y  ^  for  0  £  t  <  tg, 

thus  G(t)  £  e  t/,y  for  0  £  t  <  tg.  If  tg  <  00  then  for  t  >  tg 

G(t)  -  0  £  e"t/y.  If  F  is  NWUE  then  F(t)  >  0  for  all  t,  for  if 

F(t»)  *  0  for  a  finite  t-  then  lim  E(X-t|X>t)  -  0  <  y.  Thus 

t"*t0 

h^(t)  £  y  ^  for  all  t  £  0  and  G(t)  >  e  for  all  t  _>  0. 

The  following  inequality  (Lemma  2.2)  is  quite  an  important  tool  in 
deriving  our  subsequent  results.  It  relies  heavily  on  a  smoothing  result 
of  Feller  [3]  (Lemma  1.  p.  510). 
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Lemma  2.2.  Let  F^,  F2  be  probability  distributions  on  [0,°°)  with 
finite  means  y^  and  Assume  that  is  either  stochastically 

larger  or  smaller  than  F2,  and  that  F2  is  differentiable  with 
F2(x)  £  y^1  for  all  x  >  0.  Then: 


sup|F1(x)-F2(x) i  £  A[|yi-y2|/y1l‘ 


where  A  *  4i/6/7T. 


Proof .  By  Feller  [3],  Lemma  1  p.  510, 


(2.3) 


sup!F1(x)-F2(x)|  <2  sup|TA(t)i+^ 


where  A(x)  *»  F1(x)-F2(x) 


TA(t)  “  J  ^(t-x)VT(x)dx 


„  ,  x  1-cos  Tx 
V_(x)  -  - = - . 

ttTx 


Now,  assume  that  F^  is  stochastically  larger  than  F2»  Then: 


.00 

|  [F.(t-x)-F2(t-x)]  - 

_CO 


-cos  Tx 


[F.  (t-x)-F„(t-x)  ]  - -  dx 
»  ^TTx 


t  —  — 

-  2rr  [F,(t-x)-F2(t-x)]dx 

a.00 


mU  <VU2>' 
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Thus  from  (2.3): 


sup|f1(x)-f2(x)  |  <£  [t(p1-u2)  +-—]  . 

24 

Define  L(T)  *  T(y,-y.,)  +  - — =■,  then  a  routine  differentiation  argument 
x  ^  y^i 

gives : 

min  L(T)  -  L[ (24/Ty. (y,-y. ))1/2]  -  4/6  [l-(y_/y. ) ]1/2 
T>0  1  ^  1 

and  the  result  is  proved. 

If  F2  is  stochastically  larger  than  F^  the  analogous  result  follows 
by  similar  argument. 


3.  NBUE  Results.  Assume  that  F  is  NBUE.  Recall  that  G(t)  <  F(t) 

and  G(t)  <_  e_t/^  for  all  t  ^  0,  where  G  is  the  stationary  renewal 
distribution  corresponding  to  F.  Note  that  G* (x)  *  F(x)/y  _<  y  ^  for 
all  x.  Applying  lemma  2.2  with  F^  »  F,  F2  =  G  we  obtain: 


(3.1) 


sup|  F(x)-G(x)  1  <  (l-~)  *  P 1/2 
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(3.3)  sup|F(x)-e  x^|  £  max(sup|F(x)-G(x)|,sup|G(x)-e  X/^J)  £  p*^2  . 

Next,  by  simple  computation: 

(3.4)  sup  |  j  £  l-(yG/y)  =  p  . 

Moreover,  e  £  max(G(t)  ,e~C/^G) ,  thus: 

(3.5)  sup|G(t)-e~t/,^JG|  £  max(sup  |G(t)-e  t^|,sup|e  t/,,J-e_c/,1JG  | )  £p  , 

We  summarize  these  results  in  Theorem  3.6. 

Theorem  3.6.  Let  F  be  NBUE.  Then: 

sup|?(x)-e"x/W|  <  Ap1/2 
sup | F(x)-G(x) |  £  Ap1^2 

supjG(x)-e  X/^|  £  sup|G(B)-  y  *e  C/^dt|  £  p 

B 

sup |G(x)-e~x^Gj  £  p 

lx  l/A  O 

where  A  *  — ~  and  P  ■  1—  (V^2  ^  ^  • 

Corollary  (3.7)  below  presents  a  limit  theorem  for  NBUE  distributions. 

Corollary  3.7.  Let  {Xn,n£l}  be  a  sequence  of  NBUE  random  variables 
with  U„-EX„,  U2>„  -  EX*  and  0n  -  l-(u2n/2UJ) .  Then  Xn/un 

converges  in  distribution  to  an  exponential  distribution  if  and  only  if 


lim  pq  *  0,  in  which  case  the  mean  of  the  limiting  exponential  distribution 
equals  1. 


Proof.  The  sufficiency  of  the  condition  lim  o  *  0  follows  from 

n 

Theorem  3.6.  To  prove  necessity  assume  that  lim  Pr(X  >  ty  )  =  e  for 

n  n 

all  t  >_  0,  and  some  c  >  0.  Let  Gn  denote  the  stationary  renewal 

distribution  corresponding  to  X^,  and  H  the  stationary  renewal 

distribution  corresponding  to  X  /u  .  Then  H  (t)  *  G  (tp  )  and 

n  n  n  n  n 

»„  ■  1  -  f0  V'>dt- 

Now: 


lim  H  (t)  *  1-lim  H  (t)  «  1-lim  Pr(X  > sy  )ds  -  l-[ (l-e"Ct)/c]  . 

i CL  . Q  .  n  n  j  «  IX  IX 

n-K»  n-»O0  n-«0  J  0 


Since 

that: 


X  is  NBUE,  so  is  X  /p  ,  and  it  thus  follows  from  Lemma  2.1 
n  n  n 


H  (t)  <  e  C  for  all  n,  t  >  0 

n 


Thus  by  the  dominated  convergence  theorem: 


lim  p 

n-*» 


1-lim 

n-H»  <<0 


I  H  (t)dt  -  1-  f  [l-{(l-e"ct)/c}]dt 

Jn  n 


J0 


00  for  c  <  1 

0  for  c  *  1 

-00  for  c  >  1 


But  NBUE 


distributions  satisfy  0  <  p  <  1/2  (since  ^  >  y  by 


Chebychev's  inequality).  Thus  c  equals  1  and  lim  *  0. 
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4_.  Potential  Improvement  of  NBUE  Bound.  In  the  following  example  we  have 

a  sequence  of  IFEA  (and  thus  NBUE)  distributions  {F  ,  n  >  l),  with 

n  — 

sup  If  (x)-e  x/,yn[ 

i  » _  n  « 


It  follows  from  this  example  and  theorem  1  that  the  best  bound  of  the  form 
ct 

cp  has  a  *  1/2  and  1  £  c  £  — jj—  .  Thus  the  maximum  potential  improvement 
1/2 

in  the  bound  Ap  is  the  lowering  of  A  to  1.  This  statement  holds  for 

the  NBUE  class  as  well  as  for  the  subclasses  NBU  and  IFBA. 

The  distribution  F  is  defined  by: 

n  3 


F  (t)  - 
n  7 


J  t  <  - 
n 


e-t  t  >  -  . 
—  n 


Then: 


-1  -n 
n  +  e 


n~2  +  2(n"1+l)e_n 


-1  -1  , 
l-[(l+2n(n+l)e“n  )/2(l+ne"n  )  ] 


supjF  (x)-e  X^n|  =  1-exp [-1/ (1+ne  n  )] 


It  follows  that: 


D  *n^  +  0(n^) 
n 
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and 


pn  *  n"2  +  ^Cn"2)  • 

Thus: 

lim  CDn/|Pn|1/2]  =  1  • 
n-*30 

5_.  NWUE  Results.  Assume  that  F  is  NWUE.  Applying  Lemma  2.2  with 
*  F  and  F^  *  G  we  obtain: 

(5.1)  sup|F(x)-G(x) |  £  Ap^2  . 

We  do  not  know  of  an  analogue  of  (3.2)  for  NWUE  distributions,  but 
applying  Lemma  2.2  with  F^x)  -  e~X/^  and  F2  =  G  we  obtain: 

(5.2)  sup|G(x)-e"x/y|  £  Ap1/2  . 

Since  G(x)  £max(F(x),e  X/^p)  it  follows  from  (5.1)  and  (5.2)  that 

(5.3)  sup|F(x)-e‘x/y|  £  Ap1/2  . 

Finally  since  e  X/,p  £  min(G(t),e  and  SUpje  t^G-e  t^|  £ 

l-(y/UG)  *  p/p+1,  we  obtain: 

(5.4)  sup|G(x)-e  x^'iG|  £  Ap1//2  . 


Corollary  (3.7)  does  not  hold  for  NWUE  distributions.  While  (5.3) 

insures  that  lim  p^  =  0  is  sufficient  for  convergence  to  an  exponential 

distribution,  lim  p  =  0  is  not  a  necessary  condition.  To  see  this 

n 

consider  the  distribution  F  with  failure  rate: 

[2  0  <  x  <  1 

h(x)  -  < 

[2X*1  x  >  1  . 

Clearly  F  is  DFR,  with  finite  mean,  and  infinite  second  moment. 

Now,  for  nal,2,...  define: 

Fn(x)  -  n*1F(x) + (l-n"1)(l-e'x)  . 

Since  F  is  a  mixture  of  DFR  distributions,  F  is  DFR  and  thus 
n  n 

NWUE.  Clearly,  Fr  converges  to  an  exponential  distribution  with  mean  1. 

However  since  F  has  infinite  second  moment,  so  does  Fn,  and  thus 

p  *  00  for  all  n. 

n 

6^.  Geome trie  Sums .  Y  is  defined  to  be  a  geometric  sum  of  X  with 
parameter  p  if  Y  can  be  represented  as  X.^  with  {X^,i>^l} 
i.i.d.  as  X,  N  geometrically  distributed  with  parameter  p,  and 

N  and  {X^}  independent. 

Lemma  6.1.  If  Y  is  a  geometric  sum  of  X  with  parameter  p  then 

PY  -  ppx  where  PY  -  | (EY2/2(EY)2)-l|  and  Px  *  | (EX2/2(EX)2)-1 | . 


Proof.  Define  U  -  EX,  y2  =  EX2,  02  *  VarX,  Uy  =  EY,  U2  Y  =  EY2 
2 

and  *  VarY.  Note  that: 

(6.2)  =  li/p 

2  2-2  2-1 

(6.3)  aY  *  quS  *  +  Op  1  . 

Thus: 

(6.4)  y2,Y  ’  +  . 

From  (6.2)  and  (6.4): 

(6.5)  u2,Y/2yY  *  <I+<P2P/2U2)  * 

2  2 

We  see  from  (6.5)  that  y2  Y^2UY  —  and  only  ^  U2/2y  —  1 • 

2 

Assume  that  y2/2y  <.  1.  Then: 

(6.6)  PY  ■  i-(y2  y/2uy}  *  pU-(u2/2y2)]  *  PPX  • 

Finally  for  y2/2y2  <_  1: 

(6.7)  PY  «  (y2>Y/2y2)-l)  =  p[(y2/2y2)-l]  -  ppx 

and  the  result  is  proved. 
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Lemma  6.8.  Suppose  that  Y  is  a  geometric  sum  of  X  where  X  is  NBUE. 
Then  Y  is  NBUE.  The  analogous  result  holds  for  X  NWUE. 

Proof.  Consider  a  renewal  process  with  inter arrival  time  distribution  X. 
Binomial  sampling  with  probability  p  of  the  renewal  epochs  leads  to  an 
embedded  renewal  process  with  interarrival  time  distribution  Y.  It 
follows  that  Y  can  be  represented  as  X  +  X^,  where  X  (Y  )  has 

the  stationary  distribution  of  X(Y).  (X^}  is  i.i.d.  as  X,  N  is  geometric 

*  r  i 

with  parameter  p,  and  X  ,  tX^},  and  N  are  independent.  But,  Y  is 

N— 1 

representable  as  X  +  X.^  with  X,  N  and  {X.}  independent.  Now 

* 

if  X  is  NBUE  then  X  is  stochastically  greater  than  X  so 

Y  ■  X  +  is  stochastically  greater  than  Y  =  X  +  *  X^,  and 

Y  is  thus  NBUE.  The  analogous  argument  obviously  works  for  X  NWUE. 

Theorem  6.9.  Let  X  be  either  NBUE  or  NWUE  with  finite  second  moment. 
Suppose  that  Y  is  a  geometric  sum  of  X  with  parameter  p.  Then: 

sup jPr(Y > t)-e  tp^  |  £  A(pp)1^2 

where  A  *  and  p  *  |  (EX2/2(EX)2)-1 1 . 

Proof.  The  result  follows  from  Theorem  3.6,  Lemma  6.7  and  Lemma  6.8. 

Note  that  Theorem  6.9  applies  to  defective  renewal  processes,  which  are 
discussed  in  Feller  [3],  chapter  XI,  sections  6  and  7. 
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